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General relativistic corrections to the galaxy power spectrum appearing at the horizon scale, if 
neglected, may induce biases on the measured values of the cosmological parameters. In this paper, 
we study the impact of general relativistic effects on non standard cosmologies such as scenarios 
with a time dependent dark energy equation of state, with a coupling between the dark energy 
and the dark matter fluids or with non-Gaussianities. We then explore whether general relativistic 
corrections affect future constraints on cosmological parameters in the case of a constant dark 
energy equation of state and of non-Gaussianities. We find that relativistic corrections on the 
power spectrum are not expected to affect the foreseen errors on the cosmological parameters nor 
to induce large biases on them. 

PACS numbers: 



I. INTRODUCTION 

The complete general relativistic description of the observed matter power spectrum is, at large scales, 
significantly different from the standard Newtonian one. The observed rcdshift and position of galaxies are 
affected by matter fluctuations and gravity waves between the source and the observer, see e.g. Ref. pj. 
In addition, the matter density perturbation, 5 m , is gauge dependent while observable quantities, such as 
the power spectrum, should be gauge invariant. The standard picture looks, therefore, incomplete, and 
a general relativistic description is needed in order to correctly compute the measured observables (2, Q . 
Current observations, based on available galaxy surveys, are not affected, in practice, by general relativistic 
corrections since they appear only at very large scales. In future galaxy surveys, however, these corrections 
may interfere with the measure of other physical effects which modify the large-scale shape of the power 
spectrum. 

In this paper, we study the general relativistic effects in several cosmological scenarios, like: i) a constant 
dark energy equation of state; ii) a time varying equation of state w(a); hi) non-Gaussianities; iv) a coupling 
between dark energy and dark matter; and finally iv) massive neutrinos. For the scenarios i) and iii), 
we compute the expected errors and biases from a future Euclid-like galaxy survey by means of a Fisher 
matrix analysis, comparing the results with and without general relativistic corrections in the matter power 
spectrum. 

The structure of the paper is the following. Section |TT] summarizes the general relativistic corrections 
treatment. In Sec. IIIII the impact of general relativistic corrections in the cosmological scenarios quoted 
above is presented. The expected errors and biases on the cosmological parameters are computed in Sec. IIV 
for two particular scenarios. Finally in Sec. [V] we conclude. 



II. PRELIMINARIES 

Following the results of Refs. [H-IH, we briefly summarize the treatment of the observed galaxy power 
spectrum in rcdshift space. In linear perturbation theory, the observed (matter) density p m at a given 
rcdshift is defined as a function of the density fluctuation S m and the background (matter) density p m 

Pm EE p m (z)(l + S m ). (1) 
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For the standard ACDM cosmology, that we take as reference, the background matter density in terms of 
the current Hubble parameter Hq and today's matter density relative to the critical density Cl m reads : 

Mf) = |§an(i + 5) 3 . (2) 

In Eqs. (JT|) and (|2|), the dependence of the background density p m on the background redshift z has been 
made explicit. The ratio between the emitter and the observed frequencies at the background level is defined 
as 

v [K^u^) o a 

where and are the background photon wave vector and the background emitter/observer (e/o) four- 
velocity, respectively. At linear order in perturbation theory, the observed redshift of a given source, z, 
differs from the background one due to the matter/ gravity fluctuations that the photon encounters between 
the emitter and the observer positions. The perturbed four-velocity and photon null-vector read: 

u" = ±(l-A,v*) ; (4) 
a 



v 



8u 



K" = - \1 + —,n l + Sri 1 j , (5) 

where p = 0, ..,3 and i = 1,..,3. v l is the peculiar velocity of the observer/emitter and 8u and 8n % are 
the perturbed photon frequency and propagation direction, respectively. The conventions used for the 
perturbed Friedmann-Robertson- Walker (FRW) metric, together with a list of useful relations, can be found 
in Appendix [X] The observed (perturbed) redshift z thus reads: 



i+* s i+ 2+ fa=£i^ = (i 



Sis 

1 + I — + A+ (Bi - v l )n l 



v 



(6) 



Expressing the matter density in terms of the observed redshift, instead of the unobservable background 
one, it gives: 

Pm = Pm(z) (l + 5 m -*BE.?l)= p m ( Z ) (1 + A z ) , (7) 

\ dz p m J 

with the background matter density p m {z) function of the observed redshift z. While the density contrast 
<5 m and the redshift fluctuation 5z are gauge dependent quantities, their combination A 2 is, instead, gauge 
invariant. Notice, however, that the truly observed quantity is the galaxy number density perturbation [2|, 
0, [f| corresponding to: 

SN _ N(z)-N(z) SVo\ 

&obs = -=- — — , : = + -==- , [oj 

N N(z) Vol 

where an extra contribution from the physical survey volume perturbation appears. Being the volume density 
perturbation, <5Vol, a gauge invariant quantity, A b s is automatically gauge invariant, as it should be for any 
observable quantity. In addition, one has to introduce a bias between galaxy and matter overdensities. We 
will ignore for the moment the bias issue, deferring a brief discussion of this aspect to Sec. MI Bl 

Making use of the null energy condition and the photon geodesic equation (see Appendix [A] and also 
Refs. [H-El [f| for more details) one can write A 2 in terms of gauge invariant quantities as: 

A z = A m + 3n- V + 3(* B -$ B )-3^ * dX (* B - $ B ) , (9) 



1 See Sec. IIII Cl for non standard cosmologies in which Eq. (O is not valid. 
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where $s and arc the Bardccn potentials and A m and V are the gauge invariant matter density contrast 
and peculiar velocity, which definitions can be found in Appendix [AJ The last term in Eq. ([9]) is the usual 
integrated Sachs- Wolfe effect between the observer and the emission point with d/dX = dr+tfdi. The survey 
volume perturbation <$Vol has been carefully derived in several references, see e.g. Refs. @, H|, we therefore 
omit the details of its calculation here. Neglecting the unmeasurable monopole and dipole perturbations at 
the observer position, the expression of A f, s in terms of gauge invariant quantities reads: 



A 



obs 



2 U_ 



n- V 
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n • V + * B 



dX 

dX ( * B - *i 
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2 r° i 
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s Jo 



dX- 



-A n - $b) , 



(10) 



where r s = J s dr corresponds to the comoving distance between the source and the observer and An is the 
angular Laplacian on a unit sphere. Notice that Eq. ([TU|) holds for the standard cosmology case and reduces 
to Eq. (30) of Ref. [3| once the Eulcr equation for the gauge invariant matter velocity scalar perturbation: 



r' 



(11) 



is implemented. Also, let us emphasize that we have assumed a constant comoving source number density 
and ignored the vector and tensor contributions in Eqs. ^ and (flUl) . 

For later convenience, let us express Eqs. (O and (|10[) in the Newtonian gauge. The density perturbations 
A z and A ;, s read respectively: 
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(12) 



(13) 



where k is the lensing convergence (see Eq. (|A16I0 . and <I> jv are the scalar perturbations of the metric 
in the Newtonian gauge (see Appendix [A]) and the partial derivative d r = e\di = —n % di with e\ indicating 
the source position. With 5^ and v, we refer to the matter density and peculiar velocity perturbation in 
the Newtonian gauge. 

In the standard Newtonian approximation, the galaxy number density perturbation, A st , only gets con- 
tributions from the three first terms of Eq. (|13[) , namely from the matter density perturbation, the redshift 
space distortion term and from the convergence term. We consider that, neglecting the bias between galaxy 
and matter overdensities, the associated standard Newtonian power spectrum is related to the matter power 
spectrum evaluated in the synchronous gauge 2 in the following way: 



PA at = Pi (1 + /eff/4) 



(14) 



The latter is typically used for calculating the power spectrum when relativistic contributions can be safely 
neglected (i.e. for scales much smaller than the horizon scale). In Eq. (|14l) . the index 5* refers to the 
synchronous comoving gauge, f e g is the linear growth function and fik is the cosine of the angle between the 
line of sight and the wave vector k. In standard cosmological scenarios, the growth function / e fj is given by 
dhi 6 m /dlna. 

Notice that in Eq. (|T4|) we have ignored the contribution from the convergence term. Through all this 
study contributions from projected quantities have been neglected when computing the 3-D power spectrum 



2 For a comprehensive discussion see for example Refs. 0,0]. 
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Power spectra at redshift z=0.50 for u k =1 



Power spectra at redshift z=0.50 for |i k =0 




0.01 0.1 

k [Mpc/h] 




0.01 0.1 

k [Mpc/h] 



FIG. 1: PA obs (k) (solid lines) and PA st (k) (dashed lines) for fit = 1 (left panel) and fx^ = (right panel). The dark 
energy equation of state w is assumed to be constant and has been varied between —0.9 and —0.5. The vertical line 
corresponds to the horizon scale kn for w = —0.9. 



P(k), while their contributions have been accounted for when calculating the 2-D angular power spectrum 
Ci . Also, given that we expect galaxy formation to proceed in the potential wells of dark matter halos 3 , all 
the computed power spectra {P(k) or C£) in the following will correspond to the dark matter power spectra 
which have been obtained using the available public version of CAMB [8| . 



III. COSMOLOGICAL SCENARIOS 



We explore below the impact of general relativistic corrections in several cosmological scenarios which 
include the presence of a constant dark energy equation of state and the presence of non-Gaussianities. In 
the next section, we will estimate the foreseen errors on the several cosmological parameters involved in each 
of these two cosmologies using the Fisher matrix formalism. For the sake of illustration, we also discuss the 
effect of general relativistic corrections on the observed galaxy power spectrum in the case of a time varying 
equation of state w(a), a coupling between dark energy and dark matter and massive neutrinos. 

Unless otherwise stated, the following numerical values for the cosmological parameters have been used: 
n h h 2 = 0.02267, fldmh 2 = 0.1131, h = 0.705, the scalar amplitude A s = 2.64 x 10~ 9 and the scalar spectral 
index n s = 0.96. ^b(dm) refers to the current baryon (dark matter) energy density relative to the critical 
density and h is related to the present value of the Hubble parameter Hq = lOO/i Mpc/km/s. The sound 
speed for the dark energy fluid is fixed to c 2 = 1. 



A. Dark energy 

We first consider a cosmological model including standard cold dark matter and a dark energy fluid 
characterized by a constant equation of state w. Figure Q] shows the dark matter power spectra PA ob3 (k, Hk) 
and PA Bt (k, Hk) for both the line-of-sight (/Ltfc = 1) and the transverse (yUfe = 0) modes at z = 0.5 for several 
values of w, ranging from w — —0.9 to w = —0.5. The horizon scale, kn is also shown for the w = —0.9 case. 



We know that luminous red galaxies occupy massive dark matter halos today from weak lensing measurements • 
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Power spectra at redshift z=0.50 for n k =1 



Power spectra at redshift z=0.50 for |i k =0 
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FIG. 2: PA obs (k) (solid lines) and PA st (k) (dashed lines) for fik = 1 (left panel) and //j, = (right panel) for the 
three possible w(a) cosmologies explored here at z = 0.5. The vertical lines depict the horizon scale kn for wo — — 1 
and w a — 0. 



Notice that, in Sec. lIIIDl we discuss the k position of the dip appearing at large scales in the power spectrum 
PA ob3 (k, Hk) in the transverse direction (right plot). The modifications in the shape of the power spectra 
when relativistic effects are considered barely change when the dark energy equation of state is varied. In 
addition, the new features on the power spectrum induced by the general relativity terms appear only at 
very large scales: consequently, one would not expect much improvement on the measurement of w when 
relativistic effects are included. For the same reason, the bias induced on the dark energy equation of state 
w when the data are fitted to Pa, 4 (instead of using the full description given by Pa o68 ) is expected to be 
negligible. In Sec. IIV1 we estimate the foreseen errors on several cosmological parameters within a constant 
dark energy equation of state cosmological scenario using the Fisher matrix formalism. 

We also consider a time varying equation of state with a parameterization that has been extensively 
explored in the literature [9l4l2l|: w(a) = wq + w a (l — a). We study 3 cases: i) wq = —1 and w a = 0; 
ii) wo = — 1 and w a = 1; iii) wq = — 1 and w a = —1. For our numerical calculations, we have used the 
Parameterized Post-Friedmann (PPF) prescription for the dark energy perturbations, see Refs. [l3l - fl5l ]. 

Figure [2] shows the standard Newtonian matter power spectrum and the one with relativistic corrections 
included for the three w(a) cosmologies above, for both the line of sight and the transverse modes at redshift 
z = 0.5. Notice that the dependence of the general relativistic corrections on the time varying equation of 
state w(a) is as mild as for the case of constant w. The new features due to general relativistic corrections 
in the line of sight modes barely change when different values of wq and w a are considered. Therefore, 
the extra information contained in these general relativistic terms will poorly increase the precision on the 
measurement of a time varying dark energy equation of state. 



B. Non— Gaussianity 

In this section, we take into account a non zero bias between galaxy and dark matter overdensities. 
Following the prescription of several recent studies 0, @, [l6| and considering a linear bias relation in the 
comoving synchronous gauge, the galaxy and dark matter overdensities are related by 8^ = bS^ m . If the 
primordial fluctuations are Gaussian, it is generally assumed that this bias b is scale independent. 

Deviations from Gaussian initial conditions offer a unique tool for testing the mechanism which generated 
primordial perturbations. Non-Gaussianities are commonly characterized by a single parameter, /nl- The 
local primordial Bardeen gauge-invariant potential on large scales in the matter dominated era can be written 



G 



as [17H20J 

$NG = $G + /NL ($g - ($G» » ( 15 ) 

where $g is a Gaussian random field. The non-Gaussianity parameter /nl is often considered to be a 
constant, yielding non-Gaussianitics of the local type with a bispectrum which is maximized for squeezed 
configurations [21(. The standard observables to constrain non-Gaussianities are the CMB and the Large- 
Scale Structure (LSS) of the Universe. References [22j and [23| showed that primordial non-Gaussianities 
affect the clustering of dark matter halos inducing a scale-dependent large-scale bias [24Tj30| . 

Following Refs. [22|, H3 (see also d, [H, H3 f° r recent studies with general relativity corrections), we 
consider a scale dependent bias induced by the local non-Gaussianity of the following form 

S g= bS ln wh ere b = b G + Ab , (16) 
with 6q, a constant Gaussian bias and 

"^w-wppw' (17) 

T{k) is the linear transfer function that we have taken to be equal to unity and D(a) is the growth factor 
defined as <Wi(a) / Sd m (a = !)• The linear overdensity for spherical collapse can be considered as a constant: 
S c = 1.686 |33|. 

The resulting non Gaussian halo power spectrum is shown in Fig. [3] The standard Newtonian power 
spectrum is now obtained using 

Pa s1 = PL {be + A6 + fen*) 2 , (18) 

while the general relativity-corrected power spectra is obtained expressing Eq. (|10[) in the synchronous gauge 
and replacing <5 m by the galaxy density fluctuation defined in Eq. (|T6|) . The left (right) panel of Fig. [3] shows 
the power spectra for the line of sight (transverse) modes. Note that even in the absence of general relativity 
corrections the introduction of a negative /nl induces the presence of a dip at large scales (contrarily to the 
case of positive /nl)- This can be easily understood by studying the k dependence of the factor multiplying 
Pf m in Eq. (|18[) . In the negative /nl case, once we introduce general relativity corrections, the dip at 
large scales can become shallower (deeper) in the pk = 1 (pk = 0) case. In the case of the positive /nl 
values considered here, the presence of non Gaussianities induces an increase of the Newtonian matter power 
spectrum at k < 0.01 Mpc//i. General relativity corrections may also induce an increase of the power 
spectrum but at larger scales, k < 0.001 Mpc/h. However, non-Gaussianities dominate the shape of the 
power spectrum and make the general relativity effects totally subdominant. The shape of the non-Gaussian 
power spectrum barely changes when general relativistic effects are considered, regardless of the sign of the 
non Gaussianity parameter /nl- Therefore, we do not expect an important improvement on the precision 
measurement of the different cosmological parameters nor large biases on them in a non-Gaussianity scenario 
when general relativity corrections are included, see Sec. IIVI for a quantitative analysis. 



C. Coupled and modified gravity cosmologies 

Interactions within the dark sectors, i.e. between cold dark matter and dark energy, arc still allowed by 
observations (3^ - [50| . Constraints on coupled cosmologies as well as on modified gravity models could also 
be affected by the relativistic effects on the matter power spectrum. As an illustration, we parameterize 
the dark matter-dark energy interactions at the level of the stress-energy tensor conservation equations. 
Following the notations of [471 ] . an energy momentum exchange of the following form can be introduced: 

V^ m) , = ^ and ^» T (de)u = -Q» > (19) 

with 



Q v = ^Hpdeuf^/a or Q v = ^Hp^u^ / a , 



(20) 
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Power spectra at redshift z=0.50 for (i k =1 
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Power spectra at redshift z=0.50 for |i k =0 
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FIG. 3: PA abs (k) (solid lines) and PA st {k) (dashed lines) for fj,k = 1 (left) and /Lij, = (right) for five different values 
of the parameter /nl, for a Gaussian bias &g = 2 and z = 0.5. The vertical lines depict the horizon scale ku for 
w = -0.9. 



where M^ m ' d °^ is the cold dark matter (dark energy) four velocity and £ is a dimensionless coupling, con- 
sidered negative in order to avoid early time non adiabatic instabilities [43| . In general, coupled models 
with Q v proportional to are effectively modified gravity models. Assuming a flat universe and perfect 
measurements of 51dm h 2 , flhh 2 , and of the angular diameter distance to the last scattering surface from 
Cosmic Microwave Background (CMB) observations [5l|, the amplitude of £ is degenerate with the physical 
energy density in dark matter today, Odm h 2 . Consequently, Q,^ m h 2 should be changed accordingly each time 
£ is varied, see Appendix B of Ref. [48[ for the values of Vt dm h 2 and h considered here. 

Coupled cosmologies imply some extra terms in the expression of the gauge invariant matter fluctuation 
Eq. ©. Indeed, in the case of the coupled models studied here 



dpd 



dz 



= 3 



Pdii 



l + z 



1 + Z ' 



(21) 



which directly affects the expressions for A m and A z . In the case the gauge invariant quantity defined 
in Eq. ([7]) becomes: 



3-£ 



Pdc 



n • Vdm + (*B - $b) 



d\ ( * B - $ B 



(22) 



where we have made explicit the £ dependence of the gauge invariant dark matter density perturbation: 

A dm = <5 dm + (3 - £p dc /p dm )K/n , (23) 



where 1Z is the curvature perturbation defined in Eq. (|A5[) . In the u„ e case another extra contribution 
results from the modified Euler equation. In the gauge invariant formalism for dark matter perturbations 
(see e.g. H3) 



V dm = -W dm - W B + ^H(V dc - V dm ) 

Pdm 



(24) 



Therefore, the perturbation in the number density of galaxies in the u^ e case reads 

A^=A^ m -^(V dc -V dm )-n. 

Pdm 



(25) 
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Power spectra at redshift z=0.50 for n k =1 
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Power spectra at redshift z=0.50 for |i k =0 
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FIG. 4: The left (right) panel depicts P& obs (k) and P& Bt (k) by solid and dashed lines, respectively for coupled models 
oc u^ m and /ik = 1 (fik = 0). Different values of the coupling £ are illustrated, and the redshift is z = 0.5. In all the 
models, the cosmological parameters £ld m h 2 and h have been chosen to satisfy CMB constraints. The vertical lines 
show the horizon scale for w = —0.9. 



Figure 2] depicts the resulting matter power spectra PA obB {k, fJ-k) and PA at {k, fJ-k) for coupled models with 
an interaction term proportional to for both the line of sight and transverse modes at z = 0.5 and for 
different values of the coupling £. Notice that in coupled cosmological scenarios considered here, the growth 
function appearing in the definition of P& 3t (k, ^tj.) in Eq. (|14p is given by rilnr5(j m /rilna + CPde/Pdm 0- 
Similar results arc obtained for the case in which the coupling term is proportional to u^ c . As in the case of 
the dark energy equation of state, no strong biases are expected in constraining the coupling when these new 
general relativistic terms are included in the analysis: the shape of the different curves including relativistic 
corrections barely changes when the coupling is varied. 



D. Neutrino masses 



Consider a ACDM model plus massive neutrinos of a given energy density Q, u h 2 . We would like to 
determine if the massive neutrino energy density could affect the position of the dip appearing in the dark 
matter power spectrum PA obB (k, Hk) for the transverse modes (fik = 0). In order to simplify the discussion 
let us consider the expression of A {, s in the Newtonian gauge, see Eq. (|13[) . In the approximation in which 
all projected quantities in the power spectrum computation are neglected, the dip appears for /x/. = (i.e. 
n • v = 0) when the condition 

8% m + V N - 2$^ + 1$jv + ^ + ^ W = (26) 

is satisfied. For a specific choice of redshift and of cosmology, the factor E = (2/(r s 'H) + rl/rl 2 ) does not 
depend on the wave number k. Neglecting anisotropic stress, so that >J n = $jv, and making use of the 
Einstein equations (see Ref. [EH for the prescription used here): 

k ^ N = A U ^n a (s? + 3j(l+w a )v , (27) 
fc 2 ($Ar + tW n ) = \u 2 + w a )kv a , (28) 
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the dip position in the Fourier space as a function of Q a ,5 a ,v a can be extracted 

3 



k z = 



2(5di 



3^(E-2)- ^ 



(29) 



In the previous equations the index a runs over all the relevant fluids. In principle, different scenarios with 
different f2„ will show a dip at different wave numbers. However, this k difference will vanish when the total 
matter energy density (i.e. cold dark matter plus baryons plus the neutrino contribution) is kept constant. 
Therefore, general relativity effects can not help in extracting the values of the neutrino masses. 



IV. COSMOLOGICAL PARAMETER FORECASTS AND BIASES 



In this section we explore if the measurement of the different cosmological parameters is affected by 
relativistic corrections. We present constraints from future galaxy survey measurements, making use of the 
Fisher matrix formalism. Then, we compare the cosmological parameter errors with and without general 
relativistic corrections. 



A. Methodology 

The Fisher matrix is defined as the expectation value of the second derivative of the likelihood surface 
about the maximum. As long as the posterior distribution for the parameters is well approximated by a 
multivariate Gaussian function, its elements are given by p33l - [55| 

F aP = ^Tr [C-^aC- 1 ^] , (30) 



where C = S+N is the total covariance which consists of signal S and noise N terms. The commas in Eq. (f30| 
denote derivatives with respect to the cosmological parameters within the assumed fiducial cosmology. The 
1-tr error on a given parameter p a marginalized over the other parameters is a(p a ) = a/ (F' _1 ) qq , F _1 being 
the inverse of the Fisher matrix. In order to focus on the role played by general relativity corrections, we 
have restricted the analysis to galaxy survey data, i.e. we have not included in the analysis forecasts from the 
on going Planck CMB experiment. We exploit here an enlarged version of the future Euclid galaxy survey 
experiment, with an area of 20000 de e 2 . 24 redshift slices between z = 0.15 and z = 2.55 and a mean galaxy 
density of 1.56 x 10~~ 3 , see Refs. 5(| 57| . 



Two possible fiducial cosmologies are analyzed: i) a constant w cosmology (w denotes the dark energy 
equation of state) , and ii) a constant w cosmology with the presence of primordial non Gaussianitics (char- 
acterized by the parameter /nl)- In the analysis i), the model is described by the physical baryon and cold 
dark matter densities, £l\,h? and f2dm^ 2 , the scalar spectral index, n s , h, the dimensionless amplitude of the 
primordial curvature perturbations, A s and w. In the analysis ii), which includes non-Gaussianities, the 
model is described by ilbh 2 , Od m /i 2 , h, w, the dark energy sound speed squared c 2 and the /nl parameter. 
We have therefore fixed in this case the scalar spectral index and the dimensionless amplitude of primordial 
fluctuations, expected to be measured with excellent accuracy by the CMB Planck experiment. We follow 
a conservative approach, assuming that non-Gaussianitics are constrained exclusively from the very large 
scale halo power spectrum. 

In addition to the marginalized parameter errors, the biases induced in the cosmological parameters when 
data are wrongly fitted to the standard Newtonian power spectrum, neglecting general relativity corrections, 
are also computed. The biases in the cosmological parameters read [58| 



r)O i 1 

5 Pa = (F~ 1 ) a p £ -Qpf^T- (Olbs OU) , (31) 



where the sum runs over the bins indices in i = z, k and fj,f. in the case of the 3-D power spectrum analysis, 
i.e. O = P(z,k, /ifc), and i = z and I in the case of the 2-D power spectrum analysis , i.e. O = Cg(z). F 

(k z) 

is the Fisher matrix computed with the power spectra including general relativity corrections, O ob ' s and 
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O st ' are the general relativity and standard Newtonian power spectra respectively and <Tq% is the error 
on the power spectrum with general relativity corrections. 

In the case of analysis ii), we also have determined the shifts in the parameters {fl^h 2 , Qd m h 2 , h, w, c 2 } 
that would result when mock data generated with primordial non-Gaussianities (/nl = 20 in this example) 
are fitted to a theoretical model without them. The idea is the following: if the data are fitted assuming 
a model Mi with ni parameters, but the true underlying cosmology is a model M 2 characterized by n 2 
parameters (with n. 2 > rii and the parameter space of Mi includes the model Mi as a subset), the inferred 
values of the n\ parameters will be shifted from their true values to compensate for the fact that the model 
used to fit the data is wrong. In the case illustrated here, M 2 will be the model with non-Gaussianities and 
Mi the one without non-Gaussianities, i.e. with /nl = 0. While the first ni parameters are the same for 
both models, the remaining n 2 — ni parameters in the enlarged model M 2 are accounting for the presence 
of non-Gaussianities, i.e. /nl- Assuming a Gaussian likelihood, the shifts of the remaining ni parameters 
are given by [59| : 

56' a = -(G- 1 ) aP F K 5^ a ) = l...ni,C = ni + l...n 2 , (32) 

where G represents the Fisher sub-matrix for the model Mi and F denotes the Fisher matrix for the model 
Af 2 . In the case considered in this paper, Mi is the model without primordial non-Gaussianities while 
/nl 7^ in the model M 2 so that n 2 — n\ = 1 and dtp = <S/nl = 20 . 

B. 3-D Power Spectrum 

For details regarding the calculation of the Fisher matrix for the 3-D power spectra P(k, fik) measured by 
a galaxy survey, see Ref. [60l ]. Here we perform a binning both in k and in jUjt, considering nine bins in the 
former quantity. The minimum scale k m i n is fixed to 10~ 4 h/Mpc and the maximum scale is fixed to 0.1 
/i/Mpc. 

Table U contains the I— a marginalized errors on the cosmological parameters for analysis i), with a fiducial 
cosmology with constant dark energy equation of state w = — 1. Two results are illustrated: those obtained 
with the standard Newtonian power spectrum and those obtained with general relativistic corrections in- 
cluded. Note that the errors obtained in the standard Newtonian prescription are generally 40% smaller 
than those obtained with general relativistic one, except for the w parameter in which case the tendency 
is reversed. The biases on the cosmological parameters are also presented in Tab. |TJ Note that their size 
is always smaller than the 1-cr marginalized errors and therefore these biases will barely interfere with the 
extraction of the cosmological parameters. 

Table |TI] presents the results from analysis ii), which includes non-Gaussianities with a fiducial /nl = 20. 
Recently, the authors of Ref. [32j have shown that using methods to reduce the sampling variance and shot 
noise jrUT - l63f . a full sky galaxy survey can measure general relativistic effects. We do not exploit here these 
cancellation methods, leaving these combined techniques for a future study. 

The errors on cosmological parameters resulting from the Fisher analysis are not improved including 
general relativity corrections. This fact was not unexpected, given that for the value of the /nl considered 
in this analysis the changes in the power spectrum due to general relativity corrections are almost hidden 
by the effect of non-Gaussianities, see Sec. IIII B[ Note also that the biases are always smaller than the 
corresponding l-cr marginalized errors and therefore they will have no impact on the extraction of the 
cosmological parameters. Also, we find no significant shifts in the values of the cosmological parameters in 
any of the two prescriptions when the non-Gaussianity parameter /nl is (wrongly) assumed to be zero. We 
conclude that relativistic corrections in the 3-D power spectrum will not help in constraining the cosmological 
parameters. 

Finally, we briefly comment on the dependence of the cosmological parameter errors on the maximum 
scale considered in the analysis, k max , assuming a fiducial cosmology with a constant dark energy equation 
of state w = — 1. A larger k max will imply a larger number of modes, more information from the location of 
the acoustic peaks is available and consequently the errors will be smaller. Figure [5] illustrates the size of the 
relative errors on the different cosmological parameters considered in analysis i) versus the scale k max . Going 
from k max = 0.05/i/Mpc to k max = 0.2/i/Mpc the expected errors in fi^/i 2 , Qdmh 2 , h and w are reduced by 
a factor ~ 5 while in the case of the n s parameter its error is reduced one order of magnitude. 
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Parameter P& Bt (k, [ik) 


■PA ob3 (fc,Mfc) 


Biases 


A(fi dm /i 2 ) 


0.0035 


0.0057 


7.0 10 -5 


A{VL h h 2 ) 


0.0010 


0.0016 


-8.0 10" 5 


AA B 


0.021 


0.036 


1.1 10~ 5 


Ah 


0.010 


0.017 


1.3 10~ 4 


An s 


0.012 


0.016 


-4.3 10 -3 


Aw 


0.015 


0.010 


7.7 10" 3 



TABLE I: 1-a marginalized errors from the Euclid-like survey considered here for a fiducial cosmology with a constant 
dark energy equation of state, with a fiducial value w = — 1. The third row illustrates the biases induced in the 
cosmological parameters when general relativistic corrections are (wrongly) neglected. The error on the amplitude 
of the primordial fluctuations AA B is quoted in units of 2.64 • 10~ 9 . 



Parameter P& st (k, /j,k) 


-PA obs (fc,Mfc) 


Biases 


Shifts 


A(n dm h 2 ) 


6.2 10~ 4 


6.1 10" 4 


-4.8 10 -5 


-2.4 10 -4 


A(fl h h 2 ) 


8.7 10~ 4 


9.3 10" 4 


-6.3 10~ 5 


3.4 10 -4 


Ah 


3.5 10 -3 


3.8 10~ 3 


-3.5 10 -4 


2.0 10~ 3 


Aw 


1.3 10~ 2 


2.0 10~ 3 


-3.5 10~ 3 


1.5 10~ 2 


Ac 2 


4.0 


4.3 


1.0 


-2.5 


A/nl 


3.1 


3.1 


0.7 





TABLE II: 1-a marginalized errors from the Euclid-like survey considered here for a fiducial cosmology with a 
constant dark energy equation of state, with fiducial values w = —1, /nl = 20 and c 2 = 1. The third row presents 
the biases induced in the cosmological parameters when general relativistic corrections are neglected. The shifts in 
the cosmological parameters when /nl is set to zero but the data are generated with /nl = 20 have been computed 
including general relativity corrections. Similar results are obtained using the standard Newtonian expression. 




0.05 0.06 0.07 0.08 0.09 0.1 0.2 

kmax [h/Mpc] 



FIG. 5: Illustration of the relative 1-a marginalized errors (Ap/\p\) dependence on the scale kmax using the standard 
Newtonian prescription for the cosmological parameters p = Qbh 2 ,Sldmh 2 , h,w, A s and n s for a fiducial cosmology 
with a constant dark energy equation of state with a fiducial value w = — 1. 
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C. 2-D Angular Power Spectrum 

The 2-D Ci angular power spectrum is a projection of the 3-D quantity and therefore it implies an 
integration of the 3-D power spectrum P{k) convoluted with a window function, the Bessel transform of the 
radial selection function, see Refs. [64|[65|. Therefore, the C/s are not expected to give as much information 
on the cosmological parameters as the 3-D power spectrum P(k). For the calculations presented here we 
have computed the Ci assuming no magnification bias and a constant distribution of sources with redshift. 
For details regarding the calculation of the Fisher matrix for the 2-D power spectra measured by a galaxy 
survey, see Ref. (6(| (notice that we considered l m in = 2 and l max = 400). 

Table Hill contains the l-cr marginalized errors for analysis i), a fiducial cosmology with constant dark 
energy equation of state. We show the results when the Fisher matrix formalism is applied to the 2-D angular 
power spectrum in the standard Newtonian case and in the case in which general relativistic corrections are 
included. The errors in the two prescriptions are very similar. The biases in the cosmological parameters 
are also presented, and will have very little impact in the measurement of the cosmological parameters, as 
can be noticed from their sizes. 

Table IIVI presents the analogous but for the analysis ii) with non-Gaussianities. Notice that the errors 
are exactly the same for the two prescriptions and therefore there is no improvement in the determination 
of the cosmological parameters when the general relativistic corrections are addressed in the angular power 
spectrum. The biases induced in the cosmological parameters when the data are fitted to the standard 
Newtonian power spectrum are also presented. These biases are always smaller than the corresponding l-cr 
marginalized errors and therefore, will have no impact on the extraction of the cosmological parameters. 
Also, we find no significant shifts in the values of the cosmological parameters in any of the two prescriptions 
when the non-Gaussianity parameter /nl is (wrongly) assumed to be zero. Consequently, from what regards 
the 2-D angular power spectrum, relativistic corrections will not have any impact on future measurements 
of the cosmological parameters, even if the information contained at the largest scales becomes at reach. 

Notice that, as expected, the errors on the cosmological parameters obtained exploiting the 2-D power 
spectrum are, in general, larger than in the 3-D case. In the case of analysis i), the expected errors on w differ 
by one order of magnitude. The results of the 2-D and 3-D analysis should however roughly match in the 
limit of many narrow redshift bins. We have thus carried out a new fisher matrix analysis, decreasing the size 
of the redshift bin one order of magnitude in the 2-D analysis. In the latter case, similar marginalized errors 
are obtained when exploiting 2-D or 3-D power spectrum. In the case of the non-Gaussianity parameter 
/nl the errors are three orders of magnitude larger when using the 2-D angular information. This is due to 
the fact that the 2-D Ci angular power spectrum is essentially sensitive to modes transverse to the line of 
sight, while the 3-D P{k) power spectrum benefit from extra information from the radial modes. 



Parameter CiA Bt CiA obs Biases 
A(fi dm /i 2 ) 0.0088 0.0086 0.003 



A(Q. h h 2 ) 


0.002 


0.002 


< 10~ 4 


AA S 


0.093 


0.093 


-0.03 


Ah 


0.045 


0.045 


0.003 


An s 


0.04 


0.04 


-0.02 


Aw 


0.24 


0.24 


0.01 



TABLE III: l-cr marginalized errors from the Euclid-like survey considered here for a fiducial cosmology with a 
constant dark energy equation of state with a fiducial value w = —1. The biases in the parameters are also presented. 
The error on the amplitude of the primordial fluctuations A^l s is quoted in units of 2.64 • 10 -9 . 



V. SUMMARY 

The complete general relativistic description of the observed matter power spectrum at large scales is 
significantly different than the standard Newtonian one. The observed redshift and position of galaxies are 
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Parameter 






Biases 


Shifts 


A(ft dm /i 2 ) 0.0106 


0.0106 


-0.0001 


-0.00001 


A(n b /i 2 ) 


0.0053 


0.0053 


-0.0001 


< 10" 6 


Ah 


0.0499 


0.0499 


-0.00004 


< 10~ 5 


Aw 


0.2276 


0.2264 


0.0006 


-0.0007 


Ac 2 s 


4.886 


4.771 


-0.1527 


0.025 


A/nl 


1870 


1688 


-1189 





TABLE IV: l-o" marginalized errors from the Euclid-like survey data for a fiducial cosmology with a constant dark 
energy equation of state, with fiducial values w = —1, /nl = 20, c 2 = 1. The third row presents the biases induced 
in the cosmological parameters when general relativistic corrections are neglected. The shifts in the cosmological 
parameters when /nl is set to zero but the data are generated with /nl = 20 have been computed including general 
relativity corrections. Similar results are obtained using the standard Newtonian expression. 



affected by the different matter fluctuations and by the gravity waves between the source galaxies and the 
observer, see Refs. In this paper we have studied the role of relativistic effects in the extraction of 

different cosmological parameters with the galaxy power spectrum measurements that will be available from 
future surveys. 

We have explored the impact of such corrections in several cosmological scenarios as: constant (but w ^ 1) 
dark energy equation of state, time varying w(a) = wq + w a (l — a) dark energy equation of state, coupled 
dark matter dark energy scenario, massive neutrinos and primordial non-Gaussianities. We have performed 
a Fisher matrix analysis considering data from a future Euclid-like spectroscopic galaxy survey for two 
scenarios: one with a constant dark energy equation of state, the other with non-Gaussianities. We find 
that general relativistic corrections will not interfere neither with the extraction of the standard cosmological 
parameters (as the cold dark matter and baryon densities) nor with the measurement of primordial non- 
Gaussianities. The expected marginalized errors when relativistic corrections are included in the matter or 
halo power spectra are very similar to those obtained in the standard Newtonian case. The biases induced 
in the different cosmological parameters when neglecting these relativistic effects are also negligible. We 
conclude that the measurement of the cosmological parameters will not be compromised by the presence of 
general relativistic effects, once they will be included in the analysis. 
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Appendix A: Gauge invariant formalism 

The conventions we use are from Ref. [(37j with a few exceptions. More details can be found in Ref. [47| . 
For perturbations in a flat space time, the perturbation variables can be expanded by harmonic functions 
Y^(x,k) satisfying (V x + k 2 )Y^ = 0. In the following we focus on scalar perturbations, for which we 
define: 

(Al) 

^ (5) • (A2) 



Y [S) = 








y(S) . 


J_y(S) 




k 2 l« 
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Following Ref. [67| . the FRW metric, up to first order in perturbation theory, can be written as: 

g^dx^dx" = a 2 [-(1 + 2A)dr 2 - 2B t drdx l + ( 7ij + 2H ij )dx i dx j ] , (A3) 

where 7^ is the 3D flat metric with positive signature. The perturbations A, Bi and Hij arc functions of 
time and space and are in general gauge-dependent. Expanding the independent perturbations in the Fourier 
basis, and keeping only the scalar modes, we denote: 

A -> AY {S) ■ 
Bi BY} S) ; 
Ha -> H Llij Y^ +H T Y<fK 

In the following, for the sake of simplicity, we will omit the tilde symbols in the notation. Remember that 
all these quantities are represented by the correspondent Fourier expansion and depend only on time and on 
the 3-momentum k, while the position dependence is left only in the Y basis elements. 

Using these metric perturbations, we can now define er 9 , the shear perturbation and 1Z, the curvature 
perturbation, as 

<r g = ^(H T -kB) ■ (A4) 
11 = H L + ^H T , (A5) 
which are no gauge invariant quantities. The Bardeen metric gauge invariants are defined as [68[ : 



T~L 1 

*B = A - —Og - -&g , (A6) 



<£ B = H L + iff T - ja g . (A7) 
In the same line one can define perturbations for the energy-density for a given fluid a: 

< = ; (A8) 

rr = p (l + 5 a )« + r^ , (A9) 

where v % a is the peculiar velocity perturbation of the fluid and 5 a the fluid matter density contrast. Following 
[67| one define the following gauge-invariant quantities: 

V a = v a - ^ ; (A10) 
A a = S a ~~^, (All) 

Pa ft 

where A a is the gauge invariant density contrast for the fluid a defined in the gravity rest frame. Notice 
that = v a Y l and that in Eq. © and the following, V % refers to the gauge invariant velocity perturbation 
associated to the matter component, i.e. V 1 = V£ = V m Y l . 



1. Photon wave vector: some relations 



Here we provide several relations resulting from the null energy condition K^Kn = and the geodesic 
equations K^KY^ = useful in the derivation of the expression of gauge invariant matter density perturbation 
A z defined in Eq. (|7|)- On the one hand, from the perturbed null equation K^K^ = 0, one obtains the 
following relation between the temporal and spatial null vector perturbations: 

rtHm = % + (* fl - **) - ~ - 2H T + kB^j , (A12) 
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where d/dX = d T + n l di and we have taken into account that the background null equation imposes n l rii = 1. 
We have also used the background geodesic equation giving rise to riidjU 1 = nih 1 = 0. On the other hand, 
the temporal geodesic equation K V K. V gives the following condition: 

5 (T + »-)-(*-*-)- SS (If ^0- (A13 > 



2. Newtonian gauge 



It can be useful for comparison to make a particular gauge choice. In the Newtonian gauge, a g = and 
the perturbed metric is reduced to: 

ds 2 = a 2 [-(l + 2^ N )dr 2 + (1 - 2$ N )dx i dx i \ . (A14) 

In particular, for this gauge choice metric perturbations are given by: 

= % B = A, 

$jv = -$b = -H ■ (A15) 
The convergence k, in the Newtonian gauge, reads: 

An(*Jv + *w). ( A16 ) 

where r s is the comoving distance between the source and the observer and An = cot 6 do + d 2 + 1/ sin 6 2 d 2 
is the angular laplacian on a unit sphere. 
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